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Abstract 

The dark age of the Universe is generally pointed out as the period between the 
recombination epoch {z ~ 1000) and the horizon of current observations (z ~ 
5 — 6). The arrow of time in the cosmic history describes the progression from 
simplicity to complexity, because the present Universe is clumpy and complicated 
unlike the homogeneous early Universe. Thus it is crucial to know the nature of 
the constituents, in order to understand the conditions of the formation of the first 
bound objects. 

In this paper we analyse the chemical history of this dark age through the creation 
of the primordial nuclei to the formation of the first atoms and molecules. Then we 
will describe the consequences of the molecular formation on the birth of the proto- 
objects. In this context we will mention the important works of Dennis W. Sciama 
who influenced a large number of theorists -cosmologists and astronomers- on this 
new field of research dedicated to primordial molecules. 
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1 Introduction 

At early times the Universe was filled up with an extremely dense and hot 
gas. Various important physical processes occurred as a consequence of the 
expansion. For example, Universe cooled below the binding energies of light 
elements leading to the formation of the nuclei which recombined with the 
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electrons in order to form the first atoms. After recombination, although the 
density decrease acts against molecular formation, it turns out that the tem- 
perature is small enough for this formation occurs. 

The existence of a significant abundance of molecules can be crucial on the 
dynamical evolution of collapsing objects. Because the cloud temperature in- 
creases with contraction, a cooling mechanism can be important for the struc- 
ture formation, by lowering pressure opposing gravity, i.e. allowing continued 
collapse of Jeans unstable protoclouds. This is particularly true for the first 
generation of objects. 

Interactions between primordial molecules and the cosmic microwave back- 
ground radiation (CMBR) could be important if the molecular abundances 
are sufficient (Dubrovich 1977, Melchiorri & Melchiorri 1994, Signore et al. 
1997). In particular, a resonant molecular scattering between CMBR photons 
and molecules could, on one hand, smear out primary CMBR anisotropies 
and, on the other hand, produce secondary anisotropies. 
In this paper we will describe the chemical history of the Universe from the 
formation of the first atoms until the formation of the first objects. Thus in 
Sect. 2 we will recall the nucleosynthesis then the important influence of the 
neutrinos physics in Sect. 3. The period of recombination will be described in 
Sect. 4 which will lead to the description of the chemical chemistry or forma- 
tion of primordial molecules in Sect. 5. We will discuss in Sect. 6 the influence 
of the molecules on the formation of the first objects, and on the CMBR in 
Sect. 7. The possible reionization of the Universe could be some importance 
consequences on the chemical history of the Universe, in Sect. 8 we will describe 
it, and we will give some conclusions and prospects in Sect. 9. 



2 Primordial nucleosynthesis 

The Standard Big Bang Nucleosynthesis model (SBBN) has been analysed first 
by Peebles (1966), Wagoner, Fowler & Hoyle (1967), Wagoner (1969, 1973). 
Nevertheless the Chicago Group gave strong impulse to this field, see Yang 
et al. (1979, 1984), Copi et al (1995), Schramm & Turner (1998), Schramm 
(1998) and Olive et al. (1981, 2000), and the Oxford Pole led by Sarkar (1996, 
1999). Here we only summarize the main features of the SBBN model that we 
have already presented elsewhere (Signore & Puy 1999). 
The standard model of the very early Universe is quite simple and determined 
by the three important milestones of the early Universe: the expansion gov- 
erned by General Relativity, the particle interactions governed by Standard 
Model and particle distribution governed by statistical physics. The SBBN 
depends on just one parameter which is the baryon-to-photon ratio 7] or the 
dimensionless quantity of baryon density fit,. 

In the framework of the SBBN model at times much less than one second 
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after the big bang, the Universe was a hot (T >> 10 10 K) rapidly expanding 
plasma, where most of its energy is on the form of photons at high temperature 
and relativistic particles. The weak processes such as : 

n + v < — > p + e~ , n + e + < — > p + v , and n < — > p + e~ + 17 

maintained the ratio of neutrons to protons at its thermal value, i.e. n/p ~ 1 
(where n define the neutrons, p the protons, e~ the electrons, e + the positrons 
and v the neutrinos). 

At about one second, the temperature of the universe is around 10 10 K, the 
above weak processes became ineffective. The n/p ratio froze out at about 1/6, 
leading to the collisions between neutrons and protons which form deuterium: 

n + p — > D + 7. 

The collisions with deuterium led to the formation of the helium 3 and tritium 
nuclei: 

n + D — ► 3 H , p + D — > 3 He, 
and finally to 4 He via the reactions: 

n + 3 He — > 4 He , p + 3 H — > 4 He, D + D — > A He. 

At about 100 sec. the nucleosynthesis epoch was over, thus most neutrons 
were in 4 He nuclei, and most protons remained free while smaller amounts 
of D, 3 He and 7 Li were synthetized. The low densities, the Coulomb barriers 
and stability gaps at masses 5 and 8 worked against the formation of heavier 
elements. 

The chemical composition of the Universe remains unchanged until the for- 
mation of the first stars. The yields of primordial nucleosynthesis are shown 
as a function of the baryon density in Fig.l (from Buries et al. 1999). 
It is a considerable challenge to measure the actual primordial abundances 
because of uncertainties in measuring present day abundances and because of 
uncertainties in modeling the nuclear evolution since the big bang. The mea- 
sured abundances of helium 4, deuterium, lithium 7 and helium 3 (Buries et 
al. 1999) are: 



Y p = 0.244 ±0.002 
(D/H) P —(3A ± 0.3) x 10~ 5 
( 7 Li/H) p = (1.7 ±0.15) x 10" 
( 3 He/H) p =(0.3± 1) x 10~ 5 



Fig.l summarizes the measured abundances with their uncertainties. Since 
1980 and until recently, cosmologists introduced a concordance interval for 
the baryon density where the predicted and measured abundances for all the 
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four light elements are consistent; see for instance Copi et al. (1995) who 
derived 

0.007 < tt b h 2 < 0.024. 

A change occured in 1998-1999 with the accurate determination of the primor- 
dial deuterium abundance which led to the most accurate determination of the 
baryon density: Fig.l shows the concordance intervals for each element and 
the baryon density predicted by the deuterium measurement (vertical band). 
From the latest determined primordial D abundance (Tytler et al. 2000): 

(D/H) p = (3.3 ± 0.5) x 10~ 5 (95% c.L). 

Buries et al. (2000) give the following precise determination of the baryon 
density: 

tt b h 2 = 0.0189 ± 0.0019 (95% c.L). 

Until this year, the agreement between the SBBN predictions and observed 
abundances makes the SBBN a cornerstone of the Big Bang cosmology 
But the recent cosmic microwave background measurements of BOOMERANG 
and MAXIMA (de Bernardis et al. 2000, Jaffe et al. 2000) favor the following 
value: 

il b h 2 = 0.030 ± 0.0321°;°°^ (95% c.l.) 

which, a priori, cannot be accomodated within the SBBN model (Buries et al. 
2000). Future satellite experiments (MAP & PLANCK surveyor) will give us 
a better cosmic microwave background determination of VL b h 2 . 
Meanwhile, the SBBN model seems to be in trouble... and cosmologists search 
for non-standard cosmology or non standard big bang nucleosynthesis. 



3 Neutrinos and nucleosynthesis 

In cosmology, neutrinos have a strong impact on big bang nucleosynthesis and 
therefore important consequences on the chemistry of the universe. Recently 
an extensive litterature discussed the possible generation of sterile neutrinos 
in the early universe as a result of their oscillations with active neutrinos . 
Dennis W. Sciama was not only concerned with active neutrinos and decaying 
neutrinos since a long time, but also with sterile neutrinos. In an European 
conference (Sciama 1999), he gave a review and pointed out the importance 
of this field 2 and described the more recent developments on the sterile neu- 
trinos and its relevance for both SBBN and the CMBR. 
Thus, we briefly recall how SBBN can constrain the properties of active neu- 
trinos and in particular how helium abundance can be used to count neutrinos. 

2 he argued: students of the CMBR have now to learn something new, thanks to the 
famous announcement by the SuperKamiokande team as Neutrino 98, confirming 
the existence of an atmospheric neutrino anomaly... (see Sciama 99) 
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Fig. 1. : Relative abundances, by number: A/H, except for He in mass fraction: 
Y p . Concordance intervals for each element (2a uncertainty) and baryon density 
predicted (vertical line) by D measurements are also shown (from Buries et al. 
1999, astro-ph/9903300). 

3. 1 Limits to the number of neutrinos 

In the standard model, neutrinos possess only weak interactions, and are cou- 
pled to the intermediate bosons Z°. Measurements of the bosonic decay width 
allow to conclude that the total number of different neutrino species is 

N v = 3.07 ±0.12 



5 



while the combined fit to LEP data gives the more accurate results: 

N v = 2.994 ±0.01. 

In the framework of the SBBN model the precise epoch of nucleosynthesis was 
controlled by the number of particle degrees of freedom: more neutrino species 
imply a higher energy density and a faster expansion, i.e. more neutrons and 
more 4 He production. Therefore the A He abundance can be used to constrain 
the number of neutrinos species. For instance, Steigman et al. (1977) gave the 
SBBN limit of 7 neutrino species while the laboratory limit was about 5000 ! 
At this time, the SBBN limit was a very good constraint. 
The helium curve in Fig.l was computed assuming N v = 3. Buries et al. (2000) 
get the upper limit Nv < 3.2 while Lisi et al. (1999) give N v < 4 depending 
on what observational constraints one uses. 

More generally, the limits on N v can be translated into limits on other types of 
particles (numbers or masses) that would affect the expansion rate just before 
the nucleosynthesis. 



3.2 Sterile neutrinos 

Sterile neutrinos are neutrinos which lack the standard electroweak interac- 
tions possessed by the active neutrinos: v e , and v T . In particular, they do 
not contribute to the decay width of Z°\ they would gravitate; therefore, if 
they existed in the early universe, their gravitational effects would influence 
the rate of expansion and therefore also the SBBN. 

Different studies have been recently done and showed the role played by neu- 
trino oscillations in SBBN; in particular, there are some controversies about 
the possible generation of lepton asymmetry by oscillations between active 
and sterile neutrinos and its consequences on the observed baryon asymmetry 
of the universe, on primordial abundances (especially 4 He and also D), on 
the CMBR (Hannestad & Raffelt 1999) and on the decaying neutrinos (see 
later in section 8). Look, for instance: Bell et al. (1998), Sciama (1998, 1999), 
Mohapatra & Sciama (1998), Dolgov (2000) and references therein. 



3.3 BOOMERANG, MAXIMA and the number of neutrinos 

As already said in Sect. 2, the pre-BOOMERANG and MAXIMA consensus 
for the central value for Q b h 2 deduced from SBBN is 

n b h 2 = 0.0189 
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while the central value implied by BOOMERANG and MAXIMA measure- 
ments (Jaffe et al. 2000) is: 

tt b h 2 = 0.03. 

The difference between the SBBN and the cosmic microwave background cen- 
tral values for flf,h 2 has already triggered discussions, explanations in a lot of 
papers. Among them: Hu et al. (2000), Tegmark & Zaldarriaga (2000), White 
et al. (2000), Kinney et al. (2000), Kurki-Suonio & Sihvola (2000), Kurki- 
Suonio (2000), Hannestad et al. (2000) and Melchiorri & Griffiths (2000). 
Let us only mention that an alternative solution can be lepton asymmetry; 
for instance, Esposito et al. (2000) find the best fit to the BOOMERANG and 
MAXIMA data with N u = 9 or suggest the best compromise N v = 13 for cos- 
mic microwave background data and light elements abundances; Lesgourgues 
& Peloso (2000) find a reasonable fit with Q^h 2 = 0.028 and an effective num- 
ber of neutrinos N„ = 6. Active-sterile neutrino oscillations can also explain 
the current observations of cosmic microwave background and light elements 
abundances (di Bari & Foot 2000). BOOMERANG and MAXIMA data and 
the number of neutrinos are a very active field of study. 



4 Recombination epoch 

In the hot Big Bang picture, after the nucleosynthesis period the natural 
question concerns the possibility to form atoms by recombination of primor- 
dial nucleus with free electrons, and the second question concerns the degree 
to which the recombination is inhibited by the presence of recombination ra- 
diation. 

The reason for which these questions are central were mentionned by Peebles 
(1968) and Peebles & Dicke (1968) based, at this epoch, on the primeval- 
fireball picture of Gamow (1948) for the evolution of the Universe. 
The cosmological recombination process was not instantaneous because the 
electrons, captured into different atomic energy levels, could not cascade in- 
stantaneously down to the ground state. Atoms reached the ground state either 
through the cosmological redshifting of the Lyman a line photons or by the 
2s — Is two photons process. Nevertheless the Universe expanded and cooled 
faster than recombination could completed, and small fraction of free electrons 
and protons remained. 

4-1 Recombination of Hydrogen 

The principles of calculations of the primordial recombination have been men- 
tionned initially by Shklovskii (1967), Novikov & Zel'dovich (1967) and Takeda 
& Sato (1968). Peebles (1968) was the first to present a theory in which the 
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very complicated recombination process is reduced to simpler terms 3 . The 
methodology of the calculations consist to take into account a three level 
atom with a ground state, first excited state and continuum, represented by a 
recombination coefficient. A single ordinary differential equation is derived to 
describe the ionisation fraction (see Peebles 1968). The assumptions are: 

• Hydrogen excited states are in equilibrium with the radiation 

• Stimulated de-excitation is negligible 

• Collisional processes are negligible 

The approximations are based on the rate of transfer of energy per unit volume 
between radiation and free electrons developped by Kompaneets (1957) then 
Weymann (1965, 1966) and coupled with the expansion of the Universe. 
Peebles (1968) computed the rate of recombination of the plasma from the 
rate coefficient for recombination to excited states of hydrogen tabulated by 
Boardman (1964), and showed that the residual ionization of the hydrogen (or 
fractional ionization x e ) is below 5.3 x 10~ 5 for a flat cosmological model at 
redshift below 1500. Jones & Wyse (1985) developped two important analytic 
approximations: one, valid at the redshift z > 1500 is a modification of the 
standard Saha equation, and the other, valid for 800 < z < 1500. They found 
that the fractional ionization is roughly 6.75 x 10~ 4 at the redshift z = 555 
much higher than the Peebles values. More recently Seager, Sasselov & Scott 
(1999, 2000) have presented a refined treatment of the recombination through 
a complete code which performs appproximate calculations history 4 . 

4-2 Recombination of Helium 

We have seen that 4 He is the second most abundant element in the Uni- 
verse after hydrogen. Zel'dovich, Kurt & Sunyaev (1968) and Matsuda, Sato 
& Takeda (1969) were the first to point out that the cosmological recombi- 
nation of helium differs from that of hydrogen. The conditions for the helium 
recombination are such that in both cases it occurs in accordance with the 
Saha equation (Bernshtein & Dubrovich 1977, Dubrovich & Stolyarov 1997). 
Recombination history of early Universe consist of three stages but helium 
recombination occurs in two steps: 

• at z ~ 6000 (or when the temperature became less than 16 000 K) the 

3 by means of a number of approximations valid in the conditions of interest, so 
that the physical situation can be described by a few variables, and the recombination 
equation can be integrated by a few variables, and the recombination equation can 
be integrated by hand, without having to abandon any of the essential elements of 
the problem as written by Peebles himself in 1968. 

4 see http : //www. astro .ubc . ca/people /scott/recf ast .html . 
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first electron recombines and singly ionized Hell forms from doubly ionized 
MIL 

• Then at z ~ 2700 (or temperature close to 7300 K), ifell recombines into 
a neutral state. 

• from z ~ 1500 recombination of hydrogen took place. 

More recently Seager, Sasselov & Scott (2000) have developped an improved 
recombination calculation of H, He (with Hell and Helll) that involves a line- 
by-line treatment of each atomic level. They found that Hel recombination 
is much slower than previously thought, and it is delayed until just before H 
recombines. 

4-3 Recombination of Lithium 

The ionization potential In = 5.392 eV of Lil is shorter than that of hydrogen 
(I H = 13.6 eV). Dalgarno & Lepp (1987) suggested the existence of Li + ions 
after the period of the hydrogen recombination, and showed that the radiative 
recombination is again effective at the redshift z ~ 450. This last point is 
particularly important for the lithium chemistry as we will see. 



5 Primordial chemistry 

From the pioneer works of Zwicky (1959) on the molecular hydrogen content 
of the Universe, presented at the San Francisco meeting of the Astronomical 
Society of the Pacific in June 1959, the idea of molecular formation at the 
immediate post-recombination took form. 

The studies on the primordial chemistry (or post-recombination chemistry) 
have been the source of a tremendous increase of the litterature. The first 
complete description of the hydrogen chemistry was developped by different 
japanese groups: Hirasawa, Aizu & Taketani (1969), Takeda, Sato & Matsuda 
(1969) and Matsuda, Sato & Takeda (1969) in the context of pre-galactic 
gas clouds or formation of galaxies. The importance of the deuterium chem- 
istry through the formation of the HD molecule was suggested by Lepp & 
Shull (1983) as lithium chemistry through the formation of the molecule LiH, 
completed by the works of Dalgarno & Lepp (1987) on the helium chemistry 
and the lithium chemistry. Latter a complete comprehensive chemistry of the 
lithium, at the post-recombination epoch, was presented by Stancil, Lepp & 
Shull (1996). 

The first chemical network including the primordial molecules (such as H2, 
HD and LiH) and ions was carried out by Lepp & Shull (1984), Latter & 
Black (1991), Puy et al. (1993) and more recently by Galli & Palla (1998), 
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and Stancil, Lepp & Dalgarno (1998). The development of the primordial 
chemistry studies owe to Dalgarno an important contribution concerning the 
calculations of the reaction rates (see Lepp & Shull 1983, Dalgarno & Lepp 
1987, Puy et al. 1993 and Galli & Palla 1998 and references therein). 

5.1 Hydrogen chemistry 

After the important works of the recombination period, it was plausible to 
imagine that molecular hydrogen plays a role in the early evolution of large 
clouds. Nevertheless H 2 was not considered as a component of the pre-galactic 
medium. 

In a review article dedicated to the hydrogen molecules in astronomy Field, 
Somerville & Dressier (1966), showed that the times scales for appreciable 
amounts to form in three-body reactions or in radiative association were al- 
ways greater than the Hubble age. In the cosmological context there are no 
grains which can catalyze the reaction of formation, thus any H 2 formed in the 
uniform background is dissociated by the radiation, until the density is too 
low to produce it. Saslaw 5 & Zipoy (1969) then Shchekinov & Entel (1983) 
showed the importance of the charge transfer reactions: 

H+ + H -> H 2 + H + (1) 

initiated by the radiation association 

H + H + — > H 2 + hv. (2) 

They pointed out that H 2 ion is converted to H 2 as soon as it is formed and 
the H 2 concentration never becomes large. 

Peebles & Dicke (1968), in a scenario concerning the origin of the globular star 
clusters, showed that these clusters may have originated as gravitationnally 
bound gas clouds before the galaxies formed, and suggested that some molec- 
ular hydrogen can form, mainly by way of negative hydrogen by the reactions 
(see also Shchekinov & Entel 1983) 

H + e- -> H +hu (3) 

followed by the reaction 

H~ + H -> H 2 + e~ (4) 

5 In 1967, William C. Saslaw was a Ph.D. student of Dennis W. Sciama at the 
Department of Applied Mathematics and Theoretical Physics of the University of 
Cambridge (see the family tree in Ellis et al. 1993) 
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with the reaction rates calculated by Mc Dowell (1961) from remarks on the 
quantum theory of the negative hydrogen ion emphasized by Chandrasekhar 
(1944), in which the electrons act only as catalysts. Thus two possible ways 
of formation were pointed out. 

Takeda, Sato & Matsuda (1969) were the firsts to study the evolution of molec- 
ular hydrogen abundance in the cosmological medium (i.e. post-recombination 
Universe) in contrast with the works by Saslaw & Zipoy (1967) and Peebles 
and Dicke (1968) in which they calculated the products of H 2 in dense clouds. 
They considered the two ways of H 2 formation and electron and proton as a 
kind of catalyzer. 

The photodetachment of H~ and the photodissociation of H 2 by the back- 
ground radiation field restrict the abundance of molecular hydrogen formed 
at early stages, although the photodestruction of molecular hydrogen is neg- 
ligible. The destruction is due to collisional dissociation (Lepp & Shull 1983). 

5.2 Deuterium chemistry 

Deuterium chemistry in a early Universe could play an important role in the 
sense that it could give some explanations on the controversy observations 
of fractional abundance in high redshift Lya (Songaila et al. 1994, Buries & 
Tytler 1996). Although similar processes contribute to the formation of HD: 

D + + H -> HD + + hu (5) 

HD+ + H -> H + + H 2 , (6) 

its formation also proceeds through: 

H+ + D -> D + + H (7) 
D + + H 2 -> H + + HD. (8) 

Thus the formation of HD is carried out in two steps. Nevertheless when 
the abundance of H 2 is sufficient, the second way of formation is dominant 
(Eq. 8). Stancil, Lepp & Dalgarno (1998) presented a complete review of the 
deuterium chemistry of the early Universe. 

5.3 Helium chemistry 

The presence of helium in the early Universe gives rise to a rich assembly of 
molecular processes as pointed out by Hirasawa (1969). The main molecular 
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species containing helium is HeH + . Processes that lead to the formation and 
destruction of the molecular ion HeH + in astrophysical environments have 
been the subject of a number of investigations. Dabrowski & Herzberg (1978) 
were the one of the firsts to introduce the possibility that this molecular ion 
exists in astrophysical plasmas; it is from these initial works that Roberge & 
Dalgarno (1982) presented the formation and destruction mechanisms. More 
recently Zygelman, Standi & Dalgarno (1998) calculated the enhancement of 
the rate coefficient for the radiative association to form HeH + : 

He + H + -> HeH + + v. (9) 

The i?2~ ions and H 2 produce HeH + by the fast reactions 

H+ + He -> HeH + + H, (10) 
H 2 + He + -> HeH + + H. (11) 

Nevertheless for all of these processes, once photodissociation ceases to be 
rapid, the H 2 ions react preferentially with neutral atomic hydrogen to form 
molecular hydrogen as we have seen. 

5.4 Lithium chemistry 

The lithium chemistry is initiated by the recombination of lithium, which 
occured near the redshift z ~ 450. The formation of the molecular ion LiH + 
formed by radiative association processes: 

Li + + H -> LiH + + hv, (12) 

H + + Li -> LiH + + hv (13) 

opens the way of the formation of the LiH molecules through the exchange 
reactions 

LiH + + H -> LiH + H + , (14) 

which are more rapid than the formation by radiative association of H and Li 
atoms. 

The complete description of the lithium chemistry was emphasized by Stancil, 
Lepp & Dalgarno (1996) where new quantal rate coefficients were included. 



12 



5.5 Molecular abundances 



In order to estimate the molecular abundances which depend on the reaction 
rates 6 , it is necessary to know the thermal and dynamical evolution of the 
medium. In the Einstein-de Sitter Universe the evolution of matter tempera- 
ture T m is described by the equation: 

d'J. 'yyi - . 2 Tyyi, dci . . 

— — = ax e T r (T r - T m ) — (15) 

at a at 



where a is a constant, x e the fractional ionization of the hydrogen, T r the 
temperature of the cosmological radiation temperature and a(t) the expan- 
sion parameter. The first term describes the heat transfer from radiation to 
the electrons when the second term characterizes the cooling due to the expan- 
sion of the Universe. The result of the integration of this important equation 
led to the idea that there is thermal decoupling between matter and radiation. 
Nevertheless the formation of primordial molecules such as H 2) HD, and LiH 
can involve a thermal response due to the excitation of the rotational levels 
of these molecules. The population of the rotational levels is mainly due to 
collisional excitation and de-excitation with H, H 2 and He on one hand and to 
radiative processes on the other hand: absorption from the Cosmic Microwave 
Background Radiation (CMBR) and spontaneous or induced emission. 
The importance of radiative cooling rate by molecular hydrogen were pointed 
out by Takayanagi & Nishimura (1960), which developped theoretical calcu- 
lations of the cooling process due to the rotational excitation of H 2 molecule. 
They considered collisions with hydrogen atom in the context of interstellar 
clouds. Saslaw & Zipoy (1967) used the radiative function of Takayanagi & 
Nishimura (1960) in the context of pre-galactic gas clouds. The molecules HD 
and LiH play an important thermal role, since their allowed dipole rotational 
transitions, see Lepp & Shull (1984). Puy et al. (1993) calculated precisely the 
population of the levels and provide a molecular thermal function \& mo ; due 
to H 2 , HD and LiH. A complete description of the cooling of astrophysical 
media is found for H 2 in Le Bourlot et al. (1999), and for HD in Flower et 
al. (2000). Flower (2000) presented recently a study on the role of HD in 
the thermal balance of the primordial gas and conclude that HD is about as 
important as H 2 in the thermal balance of the primordial gas. Notice that 
the estimation of the ortho-para H 2 in the primordial gives a better estima- 
tion of the H 2 cooling as Flower & Pineau des Forets (2000) emphasized very 
recently. Furthermore some chemical reactions produce a chemical thermal 
function O c hem through the enthalpy of reaction. 



6 the complete chemical network is described in Puy et al. (1993), Galli & Palla 
(1998) and Stancil et al. (1998). 



13 



In the cosmological context the thermal evolution equation becomes 

— = ax e T r (T r -T m )- — - + — + (16) 

Puy et al. (1993) took into account the transfer process between radiation and 
matter via the Thomson scattering, coupled with the molecular source term 
and the enthalpy of the reactions in order to calculate the abundances of the 
species 7 . 

After a transient growth all abundances of the molecules become almost con- 
stant. The final abundances of H 2) HD and LiH freeze out due to the ex- 
pansion. Moreover the final values of molecular abundances depend on the 
choice of the cosmological parameters. Galli & Palla (1998) summarized and 
discussed the chemistry of early Universe for different cosmological models. 
In the standard model defined by the following cosmological parameters: 

H = 67 for the Hubble constant and 7710 = 4.5 for the baryon-to-photon 
ratio in an Einstein-de Sitter Universe. In this context we have at the end of 
the recombination, the following initial abundance: 

D/H ~ 4.3 x lfr 5 , and Li/H ~ 2.2 x 1CT 10 . 

In this context we obtained the following final abundance at z = 5: 

H 2 /H ~ 1.1 x 1CT 6 , HD/H ~ 1.2 x lfr 9 and LiH/H ~ 7.2 x lfr 20 

and for the main molecular ions with i e ~3x 1(U 4 : 

H^/H ~ 1.3 x 1(T 12 , HD+/H - 2.1 x 1(T 18 , H 2 D + /H ~ 5.1 x 1(T 14 , 
HeH + /H ~ 6.2 x 1(T 13 and LiH + / H - 9.4 x 1(T 18 . 

6 Primordial objects 

The developement of structure in the universe was well advanced at high red- 
shifts. For example, quasars have been detected nearly to z ~ 5, and the 
most distant galaxies to even greater distances. Still unknown, however, is the 
epoch during which the first generation of objects was formed. Hoyle (1953) 
were one of the firsts to suggest that the importance of the fragmentation 
process of extragalactic hydrogen cloud gas clouds, with particular reference 
to the formation of galaxies out of hot hydrogen clouds. From these pioneer 

7 The numerical integration of the coupled chemical equation is an initial value 
problem for stiff differential equations (see Puy et al. 1993). 
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works Hunter (1962) developped an instability theory of a collapsing gas cloud. 
Gravitational instability and thermal instability were supposed to be the main 
processes to form condensations in a dilute gas. In the framework of the grav- 
itational instability theory, each structure started as a tiny local overdensity; 
nevertheless very little is known about the protoclouds. Fluctuations that 
survive decoupling are subject to gravitational instabilities if they are on suf- 
ficiently large scale. The general treatment of the gravitational instability in 
an expanding Universe, developped by Lifshitz (1946) and Lifshitz & Kha- 
latnikov (1963), showed that a condensation due to gravitational instability 
cannot grow so fast during the time scale of the Universe. Parker (1953) sug- 
gested that the problem of condensation can be understood as a consequence 
of instability on the thermal equilibrium of a diffuse medium. 
Since the works of Jeans (1902) the importance of thermal instability as an- 
other possible mechanism for galaxy formation were made by Hoyle (1958) 
and above all Field (1965). Kato, Nariai & Tomita (1967) examined the ther- 
mal instability of a dilute gas in expanding Universe. Sofue (1969) found that 
formation of galaxies could be initiated by thermal instability due to radiative 
cooling, and developped by Saito (1969) in the context of a uniformly rotating 
homogeneous medium. The formalism were definitively established by Kondo 
(1970) in an expanding medium. 

Doroshkevish, Zel'dovich & Novikov (1967) showed that if the cloud contracts 
adiabatically, hydrogen atoms in it are soon collisionally ionized again and 
the cloud is finally dispersed by radiation pressure. Therefore, some kind of 
cryogen is necessary to form a bound system as mentionned Takeda, Sato 
& Matsuda (1968). The problem was to find a good cryogen and so primor- 
dial molecules offered an interesting solution. Matsuda, Sato & Takeda (1969) 
realized that hydrogen molecule could be a possible cooling of pre-galactic 
gas clouds from the works on the H2 cooling in the interstellar clouds of 
Takayanagi & Nishimura (1960). Yoneyama (1972) pointed out that after hy- 
drogen molecules are formed, the cooling of the gas by these molecules plays 
an important role on the further condensation and fragmentation of the cloud. 
Since these first works, the hypothesis that galaxies formed from density per- 
turbations which collapse within the expanding background of a Friedmann 
universe appears worthy of continued investigation. Nevertheless all of these 
first models considered the collapse of non-rotating spherical clouds. Thus 
Hutchins (1976) studied the thermal effects of H 2 molecules in rotating and 
collapsing spheroidal gas clouds, and suggested an early period of formation 
of objects in the mass range of ordinary stars; which led to the model of Silk 
(1977) concerning the fragmentation of cosmic gas clouds consisting of gas of 
a predominantly primordial composition. The contraction of a such cloud is 
initially adiabatic. 

Early studies focused on the chemical evolution and cooling of primordial 
clouds by solving a chemical reaction network within highly idealized col- 
lapse models (Palla et al. 1983, Mac Low & Shull 1986, Puy & Signore 1996, 
Tegmark et al. 1997). Some hydrodynamic aspects were studied in spherical 
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symmetry and multidimensional studies by Anninos & Norman (1996) and 
more recently by Abel et al. (2000). 

A small fractional abundance of H 2 molecules is formed via the H~ way. The 
small residual ionization remaining after decoupling is sufficient to produce 
enough H 2 to radiate away the compressional energy of collapse. Further col- 
lapse proceeds almost isothermally, and considerable fragmentation occurs. 
This process is terminated only when at sufficiently high densities the H 2 
molecules are destroyed by the reaction: 

H 2 + H + H -> 3H. (17) 

Palla, Salpeter & Stahler (1983) examined precisely the three body reactions 

H + H + H^H 2 + H (18) 
H + H + H 2 -»• H 2 + H 2 , (19) 

and investigated the thermal and chemical evolution of a collapsing spherical 
cloud composed of pure hydrogen gas, and various regimes of oscillation in the 
collapse can occurs as Lahav (1986) showed. 

Lepp & Shull (1984) emphasized that the dipole rotational transitions of HD 
and LiH are particularly important at high density and low temperature. 
In this context Puy & Signore (1997) examined the evolution of primordial 
molecules in a context of gravitational collapse and showed how the abun- 
dances can be modified. The importance of HD molecules was pointed out by 
Puy & Signore (1998b) where they analysed the ratio between the molecular 
cooling due to HD and that due to H 2 , and found that the main cooling agent 
around 200 K is HD. This results was confirmed by Okumurai (2000), Uehara 
& Inutsuka (2000) and Flower et al. (2000). 

The cooling of LiH is less important. Lepp & Shull (1984) and Puy & Signore 
(1996) considered the LiH molecules formed through only the radiative asso- 
ciation. Stancil, Lepp & Dalgarno (1996), Dalgarno, Kirby & Stancil (1996) 
and Gianturco & Gori-Giorgi (1996) proposed a new lithium chemistry with 
new ways of LiH formation. In this context Bougleux & Galli (1997) then 
Puy & Signore (1998a) developped a complete lithium chemistry in a collaps- 
ing cloud. Nevertheless at the temperature higher than 200 K, H 2 and HD 
cooling are thermically dominant as Abel et al. (1997) pointed out in three- 
dimensional numerical simulations, which lead to a scenario of the formation 
of first stars in the Universe due to the fragmentation of primordial gas (see 
also Bromm, Coppi & Larson 1999). 

Thus primordial molecules play an important role in the pregalatic gas par- 
ticularly on the thermochemical stability as mentionned by Rozenzweig et al. 
(1994) and on the problem of the minimum mass which virialized in order to 
form the first cosmological objects (Tegmark et al. 1997). 
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7 Anisotropies of CMBR and primordial molecules 

The study of the CMB anisotropies (CMBA) -primary and secondary ones- 
are important tools to study the origin and the evolution of perturbations. 
Different physical processes are responsible for the production of the CMBA: 
each one associated with a typical angular scale. 

Among the different kinds of secondary CMBA, the Rees-Sciama effect (Rees 
& Sciama 1968) can produce a detectable signal only if very large masses are 
involved (M > 10 16 M©) and can give significant contributions only at large 
angular scales (degree scale) and in Q ^ 1 scenarii. The Sunyaev-Zel'dovich ef- 
fect (Sunyaev & Zel'dovich 1972) is important at cluster scales (subarcminute 
scale) . At very small scales (subarcminute scale) two kinds of secondary CMBA 
can be sizeable: thermal emission by dust in galaxies and primordial molecu- 
lar lines produced by resonant elastic scattering. Here, we only consider the 
primordial molecular lines. 



7. 1 Interaction between the CMB and primordial molecules 

From an initial idea of Zel'dovich, Dubrovich (1977) showed that resonant 
elastic scattering must be considered as the most efficient process in coupling 
matter and radiation at high redshift. They noted that the cross section for 
resonant scattering between CMB and molecules is several orders of magni- 
tude larger than that between CMB and electrons: even a modest abundance 
of primordial molecules would produce significant Thomson scattering. At this 
point, every velocity field (due to molecular motion or to cloud infall) would 
leave its imprint on CMB via Doppler shift. This technique for exploring the 
early universe has been revised by the group of Melchiorri (De Bernardis et al. 
1993, Melchiorri & Melchiorri 1994, Maoli et al. 1994, Signore et al. 1994). A 
careful analysis of various primordial molecules has led to the identification of 
LiH and LiH + as the best candidates (Dubrovich 1977, de Bernardis 1993). 
Note only that this conclusion was relevant because all these authors assumed 
a value for the LiH abundance (i.e. LiH j H ~ 3 x 1CT 8 — 3 x 1CT 9 ) which was 
adopted -at that time- by all the experts in primordial chemistry. 
An attempt to search for LiH has even been carried out with the IRAM 30m 
telescope (de Bernardis et al. 1993). Because the resonant scattering is an elas- 
tic process, it can result in a primary CMBA attenuation and in a secondary 
CMBA production, exactly as in the case of an early reionisation of the uni- 
verse. Note that this effect is strongly frequency-dependent while the effect of 
Thomson scattering in a ionized medium is not frequency dependent. Finally, 
one can show that the blurring of primary CMBA by resonant scattering lines 
of LiH can occur if LiH/H > 3 x 10~ 8 for v < 60 GHz and for intermediate 
angular scales (see Maoli et al. 1994). 
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7.2 Possible molecular signals ? 

Just as primary CMBA could be erased, secondary CMBA could be expected 
by molecular resonant scattering. In the framework of a Cold Dark Matter 
(CDM) scenario -and also in a more general case- Maoli et al. (1996) have 
calculated the intensity the line width of the expected signal during the three 
stages of the evolution of a very simple protocloud model: linear phase, turn- 
around epoch and non-linear collapse. In particular, they have shown that for 
standard observational conditions (1CT 5 < Av/v < 1CT 4 ; AI/Icbr > 10~ 4 the 
beginning of the non-linear collapse phase of a protocloud, just after its turn- 
around epoch, is the best one for a detection of the first two LiH rotational 
lines which are observable in the frequency range 30 < v < 250 GHz with 
an angular scale range 7"< 9 < 20". They concluded that the IRAM 30-m 
telescope with its 10" of angular resolution and its frequency channels, is well 
suited for this kind of observations, see also Signore et al. (1997). 

7. 3 Discussion 

The effects of LiH molecules on the CMBA -erasing of primary CMBA, cre- 
ation of secondary CMBA- strongly depend on the final LiH abundance which 
is a function of the lithium produced by primordial nucleosynthesis and of the 
efficiency of its conversion to LiH . With the quantal rate coefficient for the ra- 
diative association of lithium and hydrogen given by Stancil et al. (1996) (see 
Sect. 6) the fractional abundance of LiH is much smaller than that of the above 
CMBA studies and therefore leads to no erasing of the primary CMBA and 
no production of secondary CMBA . But let us only emphasize that, again at 
present times, the primordial lithium abundance and the percentage of lithium 
converted to LiH are both quite uncertain ! 



8 Reionization 

In the 1960's years there were several ways which led to the idea that the 
medium at high redshift could be ionized: The radiation from a very distant 
object could be scattered by free electrons (Field 1954, Gunn & Peterson 
1965); radio waves of low frequency from extragalactic sources would be ab- 
sorbed (Field 1954, Sciama 1964a, Ericson & Cronyn 1965), and the radiation 
from radio sources could suffer dispersion (Haddock and Sciama 1965). Sciama 
(1964a, 1964b) and thereafter Rees & Sciama (1966) were one of the first to 
put forward that the intergalactic hydrogen could be partially ionized at large 
scale. 
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The absence of a Gunn-Peterson (Gunn & Peterson 1965) effect in the spectra 
of high redshift quasars implies that the intergalactic medium is highly ionized 
at low redshifts (i.e. z ~ 5). These observations suggested that the intergalac- 
tic medium was reionized in the redshift interval (5 < z < 300) 8 . Thus at 
epochs corresponding to z ~ 1000 the intergalactic medium is expected to re- 
combine and remain neutral until sources of radiation and heat develop that 
are able of reionizing it. 

A number of suggestions have been made for possible alternative sources of 
intergalactic ionizing photons. Sciama (1993a) gave an exhaustive description 
of the scenario of these alternatives sources. 

The scenario through decay of dark matter particles has been proposed by 
different authors, the first was pointed by Cowsik (1977) and extented by de 
Rujula & Glashow (1980) from massive neutrinos decayed into lighter neutri- 
nos and UV photons. The hypothesis of the photoionization of the intergalatic 
medium by radiatively decaying neutrinos or photinos has been previously 
formulated by Sciama (1982a, 1982b, 1982c), which led also to a scenario de- 
scribing the ionization of HI regions in the galaxy (Sciama 1990a) or the 
ionization of Hydrogen throughout the Universe (Sciama 1990b, 1993b, 1994, 
Sethi 1997). In this context Salati & Wallet (1984) reinvestigated the recom- 
bination of neutral hydrogen taking into account light neutral fermions, stable 
or radiatively unstable. They found a limit on the fractional ionization: 
4 x 10~ 4 < x e < 2 x 10~ 3 

instead of the previous baryon-dominated universe result: 

3 x 10~ 5 <i e <3x 10~ 4 . 

Sciama (1988) from the works of Cabibbo, Farrar, Maiani (1981), Sciama & 
Melott (1982) and Sciama (1982c, 1984) proposed a scenario of photino decay 
in order to explain the observation of the ionization of Lyman-a clouds at 
large redshifts. 

Bowyer et al. (1999) observed the spectrum of the night sky with an extreme- 
ultraviolet spectrometer covering the bandpass from 350 to 1100 A on board 
the Spanish satellite MINIS AT- 1 launched in 1997 April 21. The observed is 
far less intense than that expected in the Sciama model of radiative decay of 
massive neutrinos, but they cannot exclude the earlier model of Melott (1984), 
based also on the theory of Sciama with a decay energy which is somewhat 
greater than 13 eV, and different than the 13.6 KeV of the Sciama model. 



9 Conclusions 

First, let us recall that the latest CMB results require a higher baryon density 
than allowed by SBBN theory and observations and therefore new physics -for 

8 The detection of CMBR anisotropics rules out a fully ionized intergalactic medium 
beyond z ~ 300, see Scott, Silk & White (1995) 
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example a strong lepton asymmetry- 
Then, we have seen that the chemistry could have consequences on the forma- 
tion of the proto-objects. Nevertheless many important questions still await 
answers. In primordial chemistry a lot of reaction rates are estimated only 
and particularly those involved in the formation of LiH molecules. Recently 
Dickinson & Gadea (2000) developped a fully quantal description of the ra- 
diative association in Li + + H collisions, and showed that the rate coefficient 
is at least five orders of magnitude larger than that for the classical radiative 
association Li + H. Thus the radiative cooling due to LiH could be changed 
(LiH has a strong dipole moment). 

Calculation of the rate coefficients, for rotational transitions induced in col- 
lisions between the primordial molecules and atoms, is crucial. For example 
Roueff & Zeippen (1999) presented calculations of rotational excitation of HD 
molecules by He atoms, and provide a new estimation of molecular cooling 
and heating. In this context Le Bourlot et al. (1999) and Flower et al. (2000) 
showed the importance to use the recent quantum mechanical calculations of 
cross-sections for rotational transitions, in order to estimate respectively the 
H 2 cooling and HD cooling. Thus recently the rate coefficients are accessed 
in the software library of Collaborative Computational Project 9 , the UMIST 
data-base 10 which provide a large data of reaction rates and collisional co- 
efficients. Primordial chemistry is still in a nascent stage: the development of 
computing technics and specifical software for the quantum mechanics will 
open large possibilities. 

Primordial molecules could play a role at the turn-around period, i.e. where 
the expansion of the density perturbation is maximum and begins to collapse, 
and a modification of the temperature at this turn-around point can occur 
(Puy 2000b, Puy & Signore 2001). Moreover the primordial molecules play a 
crucial role on the evolution of the first objects, through the process of frag- 
mentation. The first objects could be massive stars. Thus, as soon as stellar 
processes occur, proto-objects can lead through SN explosions to the forma- 
tion of other molecules species, such as CO, CI or HCN. They in turn are 
important sources of contamination of the medium, and thus can offer differ- 
ent ways of early pre-biotic molecular formation. Chakrabarti & Chakrabarti 
(2000) showed that a significant amount of adenine (H 5 C 5 N 5 ), a DNA base, 
may be produced during molecular cloud collapse through the chain reaction 
of HCN addition. The existence of primordial carbon or nitrogen could also 
produce primordial bio-molecules or pre-biotic molecules just after the cosmo- 
logical recombination period (see Puy 2000a). 

The researches and the perpectives on primordial molecules are very large. 
The ESA cornerstone mission as Far InfraRed and Submillimetre Telescope 
(FIRST), or the Atacama Large Millimetre Array project (or ALMA), with a 
collecting area of up 10 000 m 2 , will offer interesting perspective of detection. 



9 http://ccp7.dur/ac.uk/ 

10 http : //saturn/ma . umist .ac.uk: 8000/ 
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